Extraction of proton form factors in the timelike region from 
single-polarized e + e — ► p p events 

Andrea Biancon 
Dipartimento di Chimica e Fisica per i Materiali e per I'Ingegneria, 
via Valotti 9, 25100 Brescia, Italy, and 
Istituto Nazionale di Fisica Nucleare, Sezione di Pavia, 1-27100 Pavia, Italy 

Barbara Pasquinjj] and Marco Radic 
Dipartimento di Fisica Nucleare e Teorica, Universita di Pavia, and 
Istituto Nazionale di Fisica Nucleare, Sezione di Pavia, 1-27100 Pavia, Italy 

We have performed numerical simulations of the single-polarized e + e~ — > p ' p pro- 
cess in kinematic conditions under discussion for a possible upgrade of the existing 
DAFNE facility. By fitting the cross section and spin asymmetry angular distri- 
butions with typical Born expressions, we can study the conditions for extracting 
information on moduli and phases of the proton electromagnetic form factors in the 
timelike region, which are poorly known and whose preliminary data show puzzling 
features. We have explored also non-Born contributions by introducing a further 
component in the angular fit, which is related to two-photon exchange diagrams. 
Using a dipole parametrization, we show that these corrections can be identified if 
larger than 5% of the Born contribution; we also explore the conditions for extract- 
ing information on the phase and, consequently, on the relative weight between their 
real and imaginary parts, which are presently unknown. 

PACS numbers: 13.66.Bc, 13.40.Gp, 13.40.-f, 13.88. +e 



*Electronic address: andrea.bianconi@bs.infn.it 
^Electronic address: barbara.pasquini@pv.infn.it 
■^Electronic address: marco.radici@pv.infn.it 



2 



I. INTRODUCTION 

The electromagnetic form factors are one of the most relevant sources of information on 
the structure of hadrons and on their internal quark/gluon dynamics. A lot of data for 
nucleons have been accumulated in the spacelike region using elastic electron scattering (for 
a review, see Ref. [lj and references therein). While the traditional Rosenbluth separation 
method suggests the well known scaling of the ratio Ge/Gm of the electric to the magnetic 
Sachs form factor, new measurements on the electron-to-proton polarization transfer in 
e~p — ► e~f> scattering reveal contradicting results, with a monotonically decreasing ratio for 

n 

increasing momentum transfer — q 2 = Q 2 |2j. This has led to several theoretical works in 
order to test the reliability of the Born approximation underlying the Rosenbluth method 
(see Ref. and references therein). 

Hadron form factors are experimentally less well known in the timelike region q 2 > 0, 
where they can be explored in e + e~ annihilations or hadron-hadron collisions (for a review 
see Ref. |4|). Being complex objects, their absolute values can be extracted by combining the 
measurement of total cross sections and center-of-mass (cm.) angular distributions of the 
final products. The phases are related to the polarization of the involved hadrons (see, e.g., 
Refs. 0,0,0]), but they have not yet been measured. The available unpolarized differential 
cross sections have always been integrated over a wide angular range, such that the relative 
weight of \Gm\ and \Ge\ is still unknown. In the data analysis either hypothesis Ge = or 
\Ge\ = \Gm\ were used for all q 2 , which are both not justified a priori. 

Nevertheless, the few available results show quite interesting properties. The amplitudes 
in the timelike and spacelike regions are connected by dispersion relations Q|; consequently, 
at sufficiently large \q 2 \ the form factors should behave in a similar way in both regions. 
A fit to the existing proton \Gm\ data for q 2 < 20 GeV 2 p( suggests that surprisingly 
this asymptotic regime has not yet been reached. Moreover, the very recent data from the 
BaBar collaboration on \Ge/Gm\ 9] show that the ratio is larger than 1, contradicting the 
spacelike results with the polarization transfer method [2] and the previous timelike data 
from LEAR jl^. Also the few neutron \Gm\ data are unexpectedly larger than the proton 
ones in the corresponding q 2 range Finally, very close to the threshold q 2 = 4m 2 (with 
m the nucleon mass), the trend of data suggests that the form factors could be affected by 
interesting subthreshold resonance structures (for more details, see Ref. jl2j). 



The ongoing discussion about the upgrade of the DAFNE facility by enlarging the 
cm. energy range from the <fi mass to 2.5 GeV while keeping a luminosity of 10 32 cm _2 s _1 , 
and, in particular, by inserting a polarimeter around the interaction region Q| , would allow 
to explore the production of polarized pp, nn, AA, and EE pairs with great precision. As for 
protons, in a previous paper ^| we explored the conditions under which the ratio \Ge/Gm\ 
could be extracted at DAFNE-2 from the angular distribution of the unpolarized cross 
section at any given q 2 . We showed that with 300 000 events in the considered 3.8 < q 2 < 6.2 
GeV 2 region, the uncertainty on the ratio is about 10%. Here, we consider the case where 
the final proton is polarized normally to the reaction plane. Using the same sample of events, 
we try to determine the conditions to realistically extract the phase of Ge/Gm from the 
simulation of a single-spin asymmetry measurement. 

We also made some simulations of non-Born contributions, mainly related to two-photon 
(27) exchange diagrams, which seem to play a crucial role in the analysis of the Ge/Gm 
ratio in the spacelike region Q]. Theoretically, these mechanisms are very poorly known; 
in particular, the relative size of the real and imaginary parts is undetermined. Therefore, 
in our simulations we explore the most favourable conditions to set upper limits on their 
magnitude and verify if their contribution can be measured at DAFNE-2. In Ref. [15], we 
assumed the 27 amplitude mainly real, and we came to the preliminary conclusion that 
their effect can be isolated only if their absolute size is at least 5% of the Born contribution. 
Here, at variance, we check that if the 27 amplitude is mainly imaginary, its effect could be 
detected more easily in single-spin asymmetries than from unpolarized cross sections. 

In Sec. |Hj we briefly review the necessary general formalism. In Sec. Hill we outline the 
main features of our Monte Carlo simulation of the e + e~ — > p p process. In Sec. IIV[ we 
describe the steps to extract information on form factors by reconstructing the simulated 
normal polarization with a two-parameter fit. In Sec. |Vj we discuss the results. Finally, 
some concluding remarks are given in Sec. IVII 



II. GENERAL FORMALISM 

The scattering amplitude for the reaction e + e~ — > pp, where an electron and a positron 
with momenta k\ and k 2 , respectively, annihilate into a proton-antiproton pair with momenta 
Pi and P2, respectively, is related by crossing to the corresponding scattering amplitude for 
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elastic e~p scattering. There are several equivalent representations of it; here, we use the 
one involving the axial current following the scheme of Ref. [16]. The scattering amplitude 
can be fully parametrized in terms of three complex form factors: G e{q 2 ,t) , G a/(<? 2 ,t) , and 
Ga(q 2 , t), which are functions of q 2 = (A; 1 + /c 2 ) 2 and t = (k 2 —pi) 2 . They refer to an identified 
proton-antiproton pair in the final state, rather than to an identified isospin state. In the 
following, we will use cos 9 instead of the variable t, with 9 the angle between the momenta 
of the positron and of the recoil proton in the cm. frame. 

In the Born approximation, Ge and Gm reduce to the usual Sachs form factors and do 
not depend on cos 9, while Ga = 0. For the construction of our Monte Carlo event generator, 
we rely on the general and extensive relations of Ref. pij] up to the single polarization case. 
These relations assume small non-Born terms and include them up to order a 3 (with a the 
fine structure constant), i.e. they consider 27 exchanges only via their interference with 
the Born amplitude. These corrections introduce explicitly six new functions that all can 
depend on cos 9: the real and imaginary parts of AGe and AGm, i-e. of the 27 corrections 
to the Born magnetic and electric form factors, and the real and imaginary parts of the axial 
form factor Ga- 

The role of these corrections in the unpolarized and single-polarized cross sections is 
very simple: it can be deduced from the Born term by adding the contribution of Ga 
and substituting the Born form factors with the "27— improved" ones, i.e. Ge,m{q 2 ) — * 
G EM {q 2 ) + AG E ,M{q 2 ,cosO). In Ref. Q, this fact is somehow hidden by neglecting terms 
of order a 4 . For sake of simplicity, we keep 27 effects via the axial form factor only. Within 
this scheme, when summing events with positive and negative polarization, the unpolarized 
cross section can be written as 

rlrr° 

a{q 2 )[l + R{q 2 ) cos 2 9} - b{q 2 )Re[G M {q 2 ) G A *{q 2 , cos 9)} cos9 , (1) 



dcos9 



«(9 ) = 7T~2 \ M^l + \Ge\ 2 ) , Kq 2 ) = — 2 , (2) 



2q' 2 r V t ^ 1 m 1 1 " 1 J ' q 2 r 

= T\G M (q 2 )\ 2 -\G E (q 2 )\ 2 ±_ 
{q 1 " r\G M (q 2 )\ 2 + \G E (q 2 )\ 2 'Am 2 ' 1 ] 

When the proton is polarized, the cross section is linear in the spin, i.e. da = da° (1 + 



VA), with da° from E 
states are observable 



TJ and A the analyzing power. In the cm. frame, three polarization 
6]: the longitudinal V z , the sideways V x , and the normal V y . The 



first two ones lie in the scattering plane, while the normal points in the pi x k 2 direction, the 
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x, y, z, forming a right-handed coordinate system with the longitudinal z direction along the 

77 of the outgoi,lg prot0 ' 1 - The p " is partlcularly interesti ff rf nce 14 is the ollly 

observable that does not require a polarization in the initial state [5|, |6(. With the above 
approximations, it can be deduced by the spin asymmetry between events with positive and 
negative perpendicular polarizations: 

1 do^ — do^ 



x v 



A v do' + do- 



2 J^l d(j0 sin#{cos£Im [C M (q 2 ) G E (q 2 



<T- 1 



Im 



G E (q 2 )G A *(q\ cos 9)] } . (4) 



This spin asymmetry is produced by the mechanism pi x k 2 • S#, which is forbidden in the 
Born approximation for the spacelike elastic scattering (|. 

Measurements of the unpolarized distribution (JJJ and spin asymmetry (jl} at fixed q 2 
for different 9 allow to fit the different angular terms, from which we may extract absolute 
values and relative phases of Ge/Gm, and some information on GA{q 2 , cos9). The Born 
contributions to do° and V y have a typical cos 2 9 and sin 29 behaviours, respectively, any 
deviation due to non-Born terms. The normal V y vanishes at end points 9 = 0, it and at 
threshold q 2 = 4m 2 (r =1). Interestingly, at 9 = 7r/2 the Born contribution vanishes, and 



V y gives direct insight to the 27 amplitude 

However, we note that the measurement of V y alone does not completely determine the 
phase difference of the complex form factors. By defining with Se and 5m the phases of the 
electric and magnetic form factors, respectively, the Born contribution is proportional to 
sin(<5M — 5e), leaving the ambiguity between (5m — 5e) and n — (5m — 5e)- Only the further 
measurement of V x can solve the problem, because V x oc Re(GMG* E ) oc cos(5m — Se) il- 
But at the price of requiring a polarized electron beam. 



III. GENERAL FEATURES OF THE NUMERICAL SIMULATIONS 

We consider the e + e~ — > p p process. Most of the details of the simulation are mutuated 
from a previous work Q|. Therefore, events are generated in the usual variables q 2 , 9, (f> 
(the azimuthal angle of the proton momentum with respect to the scattering plane), and S y 
(the proton polarization normal to the scattering plane). Then, the distribution is integrated 
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upon <p but not summed upon the spin. 

The exchanged timelike q 2 is fixed by the beam energy. Only the scattering angle 9 is 
randomly distributed. For a given q 2 , the Born terms in the unpolarized cross section ((!} and 
spin asymmetry (j3J) are responsible for the [1 + R(q 2 ) cos 2 9} and sin 29 Iiu.[Gm {q 2 ) G* E (q 2 )] 
behaviours, respectively. Since with the discussed approximations no other dependence in 
9 is present, observed systematic deviations from these behaviours will be interpreted as 
a clear signature of non-Born terms. The present lack of knowledge on such mechanisms 
forces us to further simplify the picture and to neglect the cos 9 dependence of the axial form 
factor by taking the first term of its expansion in powers of cos 9 (see Ref. [l^ for further 
details). 

An overall sample of 300 000 events has been considered with 3.8 < q 2 < 6.2 GeV 2 and 
| cos#| < 0.9. We recall that the integrated cross section for e + e~ — > pp in the considered 
region is approximately 1 nb; hence, at the foreseen luminosity of 10 32 cm _2 s _1 for DAFNE- 
2 [3,01 this sample can be collected in one month with efficiency 1. The lower q 2 bin does 
not include the pp threshold, because this region is characterized by peculiar mechanisms 
like, e.g., the Coulomb focussing and possible subthreshold resonances. As it will be evident 
later, the spin asymmetry between events with positive and negative S y polarization is very 
small for the first q 2 bin and would suggest to exclude it from the analysis. However, this 
would imply to reduce the size of the sample by almost a factor 2. In fact, the event 
dfcttbutt* app— l y MU lU. A Q, >eav Ulg tbe b in s a t higher f scarce* 

populated. For more details, we refer to the discussion in Ref. [15]. Here, we just remark 
that the interesting quantities related to the spin asymmetry (£Q) increase approximately in a 
linear way from threshold to q 2 ~ 5 — 6 GeV 2 : while absolute errors increase with q 2 because 
of less populated bins, the relative errors do not (see discussion in Sec. |^) . Contrary to the 
unpolarized case, there is no reason to asymmetrically split the beam time in order to make 
the statistics of each q 2 bin more homogeneous. In any case, in the following we will use the 
"default conditions", as they were defined in Ref. |l5j: the q 2 range is divided in 6 equally 
spaced bins with width Aq 2 = 0.4 GeV 2 ; for each of them, the solid angle | cos#| < 0.9 is 
divided in 7 equally spaced bins with width A cos 9 ~ 0.257. For each of the 42 bins, events 
are further divided according to their polarization, and the spin asymmetry is constructed 
according to Eq. (£Q) by taking the ratio between the difference and the sum of events with 
positive or negative polarization. 
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Events can be generated only by inserting specific parametrizations of the proton form 



factors in 
region 



;he cross section. For Ge and Gm-, several models can be considered in the timelike 
Q Q, m „s, y derived tan —ions tan the spaced r e gi o. As in Re, Q, 
we select the parametrizations of Refs. [l9| , because they have been recently updated in 
Ref. |4| by simultaneously fitting both the spacelike and timelike available data. Moreover, 
both cases release separate parametrizations for the real and imaginary parts of Ge and Gm, 
as they are needed in Eqs. and (j3J). We indicate the former as the IJLW parametrization 
(from the initials of their authors), and the latter as the Lomon parametrization. 

Since the explored q 2 range is not large, it is reasonable to make also much simpler, but 
equally effective, choices. In the previous work [15[ , we used parametrizations where all form 
factors are real and proportional to the same dipole term l/(l+q 2 /q 2 ) 2 , with q 2 = 0.71 GeV 2 , 
as in the spacelike region. The actual parameters were reduced to the ratios r e = \Ge/Gm\ 
and r a = \Ga/Gm\- However, this choice leads to a vanishing spin asymmetry, since all form 
factors have the same phase. Here, we generalize the case aiming to alternatively exphasize 
each one of the two contributions in Eq. (0J), and analyze the conditions for their extraction. 
We define the parametrization Diplz as the one with r e = 1, r a = 0, and where the relative 
phase of Ge with respect to Gm is P e = ir/2, i.e. both Gm and Ge have a dipole trend in 
q 2 but the former is real while the latter is purely imaginary; there are no effects from 27 
exchanges because Ga = 0. With Dipli, only the Born term survives in Eq. and its effect 
is maximized. The second choice is named Dip2^i and it is defined by r e = 1, r a — 0.2, and 
P e = 0, p a = 7r/2, with P a the relative phase of Ga with respect to Gm- Namely, Ge and 
Gm are equal and real; Ga has the same dipole trend, but it is purely imaginary and with 
a modulus 5 times smaller than \Gm\- In this case, in Eq. (j3J) only the non-Born term is 
active and is maximized. 

The Dipli parametrization must not be considered realistic but rather an exercise to set 
the upper limit on ImfG^]. In fact, for the considered q 2 range it is reasonable to expect 
\Ge\ ~ since the two moduli are equal at threshold. Moreover, Im[G^] cannot be 

large considering that the absorption induced by the p — p rescatterings is much below the 
unitarity limit (see Ref. and references therein). Indeed, the Born contribution to the 
asymmetry induced by the Lomon parametrization, turns out to be not larger than 10%; 
at these q 2 , the IJLW parametrization gives even a vanishing asymmetry. On the other 
side, we notice that nothing can be presently said about the relative weight of RefG^] and 
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ImjCr^]. If for the considered q 2 range we can approximately take \m\GE/u\ ~ 0, we deduce 
from Eq. (0) that the RefG^] can be best extracted from the unpolarized reaction and it 
was analyzed in Ref. Q|. Viceversa, the assumption in Dip27« of a large Ihi[Ga] selects the 
normal polarization of Eq. (j3J) as the best observable to study 27 exchanges. In this sense, 
the choice Dip27z is also an attempt to set a possible upper limit to non-Born contributions. 
However, as in Ref. |3] we stress that we set r a = \Ga\/\Gm\ = 0.2 to account for these 
effects as corrections to the Born result. In particular, we note that a correction of 20% 
does not contradict the experimental findings on Ge/Gm m the spacelike region with the 
polarization transfer method Q|. 

Event Distribution | 
25000 
20000 
15000 
10000 
5000 

°-0.8 -0.6 -0.4 -0.2 -0 0.2 0.4 0.6 0.8 
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Figure 1: Angular distribution of 300 000 events for the e + e~ — * pp process at 3.8 < q 2 < 6.2 
GeV 2 and |cos#| < 0.9, and with the Lomon parametrization of proton form factors (see text). 
The points correspond to 7 equally spaced cos 6 bins for each of 6 equally spaced q 2 bins. Solid 
curves are the results of a 3 parameters angular fit with B tx r a = \Ga/Gm\ = (see text). The 
highest curve corresponds to the lowest bin 3.8 < q 2 < 4.2 GeV 2 ; next lower curve to the adjacent 
bin 4.2 < q 2 < 4.6 GeV 2 , and so on. 




IV. RECONSTRUCTION OF THE SPIN ASYMMETRY 



In Ref. [15j , we ana lyze d in detail the conditions for an unpolarized e + e — ► pp measure- 
ment at DAFNE-2 Q] m order to estimate the achievable precision in the reconstruction 



9 



of r e = \Ge/Gm\ and in a possible detection of 27 contributions. As it has been recalled in 
Sec. HI the complete determination of Ge and Gm in the timelike region requires the knowl- 
edge of their relative phase, which in turn involves the polarized e + e _ measurement. 

Our analysis consists of five steps. For sake of simplicity, we describe it in the following 
for the case of the Lomon parametrization: 

1) We generate 300 000 events for 3.8 < q 2 < 6.2 GeV 2 and |cosfl| < 0.9 (i.e., for 
25° < 9 < 155°). The sorted events are summed upon the positive and negative 
polarizations S y and are divided into 6 equally spaced q 2 bins of width Aq 2 = 0.4 
GeV 2 . In each q 2 bin, the events are divided into 7 equally spaced cos 6 bins with 
width Acos# ~ 0.257. In Fig. Q the points show the 6 unpolarized distributions 
in q 2 , each one consisting of 7 points describing the distribution in cos#. From the 
discussion in the previous section, the highest cos# distribution corresponds to the 
lowest bin 3.8 < q 2 < 4.2 GeV 2 ; the next lower distribution to the adjacent bin 
4.2 < q 2 < 4.6 GeV 2 , and so on. 

2) For each q 2 bin, the cos 6* dependence of the unpolarized event distribution is fitted by 
a function of the form 

C q (cos6) = AN q (cos6) = A(l + R cos 2 6 - B cosfl) , (5) 

where the three fitting parameters A, R, and B, change for each different q 2 bin. The 
results of each fit are represented by the solid curves in Fig. [0 The parameter A is 
related to the total number of collected events and it does not play any role in the 
following analysis; R is connected to the angular coefficient R{q 2 ) of Eq. (0J), and B 
to the correction from 27 exchanges (for a thorough discussion, see Ref. (15j). 

3) For each q 2 , 6 bin, the events generated at step 1) are divided in two groups corre- 
sponding to positive (U) and negative (D) normal polarizations S y . Then, the spin 
asymmetry is manually constructed by taking the ratio (U — D)/(U + D). In Fig.0 the 
full squares show the cos 9 distribution for the lowest q 2 bin, while the empty squares 
refer to the largest q 2 bin. 

4) After fixing N q (cos8), for each q 2 bin the 9 dependence of the spin asymmetry is fitted 
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Figure 2: Normal polarization for the e + e 



pp process in the same conditions as in Fig. ^ Full 



squares refer to the lowest bin 3.8 < q 2 < 4.2 GeV 2 , empty squares to the highest 5.8 < q 2 < 6.2 
GeV 2 . Solid lines are the result of a 2 parameter angular fit (see text). 

by a function of the form 



where the fit parameters A% and A 2 refer, respectively, to the Born term Im[G A/ G^] 
and to the non-Born term ImfG^G^] in Eq. (J1J). The results of the fits for the lowest 
and largest q 2 bins are represented by the solid lines in Fig. 

5) The procedure is repeated 10 times, starting from statistically independent samples, 
and average values and statistical variances for A\ and A 2 are built for each q 2 bin 
and displayed in Fig. El We checked that stability of the results is ensured already 
after 6 repetitions. In the figure, each point is accompanied also by two other error 
bars, whose meaning is the following. For each repetition, A\ and A 2 are determined 
with some error from the fitting procedure. From left to right, the second error bar 
associated to each point is the average of the fitting errors upon the 10 repetitions. 
We recall that this error gives the integral deviation of the test function from the 
data points; therefore, it must be considered as a systematic error and, in principle, it 
has nothing to do with the statistical one. Finally, the rightmost error bar associated 
to each point gives the uncertainty on A\ and A 2 deriving from the uncertainty on 
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Figure 3: Values of the fitting parameters A\ and A 2 extracted from the normal polarization of 
Fig. [21 according to Eq. ©. For each point, the leftmost error bar gives the statistical error, the 
central error bar gives the fitting (systematic) error, the rightmost error bar is the uncertainty 
associated with the fitting error on the unpolarized event distribution. 

the separately fitted parameters R and B: when the latter ones oscillate between the 
extreme values allowed by their fitting error bars, the averages values of A\ and A2 
can oscillate inside the rightmost error bars displayed in Fig. El 

The Lomon parametrization includes only Born contributions to the complex Ge and 
Gm, he. Ga = in Eqs. ((H) and (JH). Consequently, the unpolarized angular distributions 
of Fig. ^ are symmetric functions in cos 9 and the spin asymmetry of Fig. 121 displays the 
typical sin 29 trend, as it should be. Consistently, the fit parameter A 2 , related to the 
interference between Ga and Ge, turns out to vanish in Fig. El within the statistical (and 
fitting) uncertainty. 

The fits of the unpolarized distribution and of the spin asymmetry are separately per- 
formed, because each single generated event can be collected in the group with positive {U) 
or negative (D) normal polarization, denning the asymmetry as (U — D)/(U + D). Our 
procedure corresponds to separately fit the denominator U + D (with the 3 fit parameters 
A,R, and B) and the numerator U — D (with the 2 fit parameters A\ and A 2 ). The un- 
certainty in the former fitting procedure does not affect much the latter one. This can be 
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visualized by the rightmost error bar associated with each point in Fig. El We checked that 
this is not accidental by repeating the fit of the spin asymmetry but artificially assuming 
the denominator U + D independent from cos 6: the results on Ai and A 2 stay basically 
the same. It must be noted, however, that in the considered q 2 range both the angular 
coefficient R(q 2 ) of Eq. Q and the non-Born contribution proportional to Im[G M G^] are 
small, because \Ge\ ~ \Gm\ an d Ga is assumed to vanish in the Lomon parametrization. In 
general, for larger q 2 the first condition may not be satisfied, while the second one appears 
too drastic. Consequently, the analysis of error bars would become more involved; but a 
discussion along this line is beyond the scope of this paper. 

As previously described, the first and second leftmost error bars in Fig. |3] represent 
the statistical and fitting errors on A\ and A 2 . The fitting error is usually a systematic 
error, except in those situations where in the limit of a large number of events the test 
function reproduces exactly the event distribution. This may happen when the number of 
fit parameters is equal to the number of bins (which is not our case, because we have 2 
parameters against 7 bins in cos#), or when the test function is the correct guess for the 
event distribution. In general, this is not possible since, in reality, the 9 dependence of the 
form factors is unknown. But this is the case of our simulation, where we have neglected 
the 9 dependence in all form factors and the fitting functions © and Q reproduce the 
complete angular dependence of the event distributions. Therefore, in the limit of a large 
number of events the statistical and fitting errors become equal; for the 300 000 events 
considered here, they are very similar, as it is evident from Fig. EJ The main consequence 
is that, in the conditions previously discussed, the statistical and fitting errors are two 
different ways to estimate the same error, and they must not be summed as one usually 
does for statistical and systematic errors. This statement can be safely generalized when 
the analysis is performed in the Born approximation, because Ge and Gm depend only on 
q 2 and the fitting functions (j3J) and © with B = and A 2 = 0, respectively, are certainly 
correct. It is more questionable when non-Born contributions are included, either in Ge/m 
or via Ga, and the cos 9 dependence should be considered also inside the form factors. 
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Figure 4: The same as in Fig. |HJ but for the IJLW parametrization (the rightmost error bar for 
each point is not displayed). 

V. DISCUSSION OF RESULTS 



In the following, we discuss the output of repeating the procedure outlined in the previous 
section but for different parametrizations, as they have been described in Sec. IIHI For sake 
of brevity, we will omit the figures with the angular distributions and will show only the 
dependence on q 2 of the parameters A\ and A 2 extracted from the fit. For each point, only 
the statistical and fitting errors will be displayed, since the sensitivity to the uncertainty on 
the R and B parameters of Eq. (JHJ is small (see Fig- ISJ) • 

In Fig. HI the average values of A% and A 2 with statistical and fitting errors are displayed 
for the considered q 2 bins for the IJLW parametrization. The A 2 is consistently compatible 
with zero because no 27 effects are included. Surprisingly, the A% is also compatible with 
zero, while at much larger q 2 this choice gives very large asymmetries, similarly to the 
Lomon parametrization jj]. Luckily, the relatively small error bars allow to unambigously 
distinguish the two realistic parametrizations, except maybe for the lowest q 2 bin. 

It is not accidental that the cleanest comparison is possible for the largest q 2 here consid- 
ered. At the threshold q 2 h = 4m 2 , we have Ge = Gm and the Born term in Eq. (0J) vanishes. 
For q 2 not much larger than q 2 h , we may assume that 



Ge ~ G 



M 



1 - a (q 2 - q 2 th ) + 



A? (e 2 -<4)+- 



(7) 
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The Born term of the spin asymmetry will behave like — sm[/3(q 2 — q 2 h )] [1 — a(q 2 — q 2 h )] w 
— P(q 2 — q 2 h ) f° r small q 2 — q 2 h . Hence, it is easy to justify the approximate linear dependence 
for q 2 < 5 GeV 2 displayed both in Fig. |3] and Fig. HJ 



ED 
















-0.4 


















-0.45 








ll « 


T 1 


T 
I 






-0.5 




si 






1 ! 








-0.55 


















-0.6 






















4 


4.5 


5 


5.5 


6 






















0.05 























♦I 


-M — 




b+ 


T 
















1 






-0.05 






















4 


4.5 


5 


5.5 


6 


q 2 



Figure 5: The same as in Fig. |21 but for the Dipli parametrization (the rightmost error bar for 
each point is not displayed). 

As anticipated in Sec. IIII[ we have also considered the Dipli parametrization which 
emphasizes the Born contribution to V y by taking Ga = and Ge purely imaginary with 
the same modulus as Gm- The corresponding fit parameters A\ and A 2 are displayed in 
Fig. 03 A 2 is obviously consistent with zero, while A\ is roughly three times larger in size than 
in the more realistic parametrizations Lomon and ILJW. Evidently, it must be considered 
as an upper limit. 

Similarly, the parametrization Dip2'ji emphasizes the non-Born contribution in V y by 
taking Ge and Gm real and Ga purely imaginary. As already discussed in Sec. I11H this is 
not an upper limit for the 27 effects, since we have imposed |GU|/]Gm| = 0.2, consistently 
with the findings in the spacelike region about the Ge/Gm ratio [2]. But it represents the 
most favourable situation to explore ImfGU] using the spin asymmetry. In Fig. 13 A\ is 
obviously consistent with zero, and A 2 indicates a clearly detectable effect. 

Analogously to the Born case, the 27 effects seem better visible at the largest q 2 here 
explored, even if the latter correspond to the least populated bins. We note also that while 
the Born term in V y approximately falls like \Ge/Gm\ for increasing q 2 and, therefore, is not 
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Figure 6: The same as in Fig. |31 but for the Dip27z parametrization (the rightmost error bar for 
each point is not displayed). 

n 

expected to fall fastly unless far beyond 10 GeV 2 jSj, viceversa, from Eq. Q together with 
Eqs. fll|2jl . the weight of the 27 contribution should fall at least as 1/y/r ~ 1/q. Therefore, 
unless the 27 effect increases with q 2 , the available range of its measurability is reduced to 
q 2 < 10 GeV 2 . 

VI. CONCLUSIONS 

We have performed numerical simulations of the single-polarized e + e~ —>■ p p process 
using a sample of 300 000 events in the kinematical region with 3.8 < q 2 < 6.2 GeV 2 , 
distributed over 6 equally spaced bins. For each q 2 bin, event have been further distributed 
over 7 equally spaced bins in cos(6'), with | cos 6*| < 0.9. For each q 2 ,9 bin, the events 
have been separated in two groups according to their positive (U) or negative (D) normal 
polarization of the recoil proton with respect to the reaction plane. 

The angular distributions of U + D (the unpolarized cross section) has been fitted with a 
three-parameter function of the type A [1+R cos 2 9 — B cos 9], where from the parameters R 
and B information on r e = \Ge/Gm\ and r a = \Ga/Gm\ can be deduced, respectively. The 
angular distribution U — D (proportional to the spin asymmetry) has been fitted with the 
two-parameter function sin 9 [Ai cos 9 — A 2 ] , where the two free parameters A\ and A 2 are 



16 



related to lm[G M G* E ] and lm[G E G* A ], respectively, i.e. to the relative phases of the complex 
form factors. 

Using realistic parametrizations for Ge and Gm deduced by fitting the available spacelike 
and timelike data, the statistical uncertainty makes it possible to unambiguously distinguish 
among different models, which are equivalent in the spacelike domain. Therefore, both 
moduli and phases can be extracted with a two-step fitting procedure, where the fitting 
errors of the first step do not affect much the uncertainty in the second step. Surprisingly, 
the extraction of phases seems better performed for the higher q 2 bins, because the spin 
asymmetry is small for q 2 close to the threshold and the relative error becomes smaller and 
smaller for increasing q 2 . 

The additional contribution related to two-photon exchange diagrams shows up as a 
deviation from the expected angular trend of the Born contribution, via the parameters 
B and A 2 . In our analysis, we modeled such contribution in terms of the complex axial 
form factor Ga, constraining its modulus not to contradict the findings in the spacelike 
region. Nothing can be said about the relative weight of RefG^] and ImfC^]. However, 
assuming that for the considered q 2 range Ge and Gm 

are similar and with small imaginary 
parts, information on RefG^] and Im[GU] can be extracted from the non-Born term of the 
unpolarized cross section and of the spin asymmetry, respectively. Using a simple dipole 
parametrization, our finding is that the non-Born term can be identified if its size is at least 
5% of the Born contribution; if ImfG^] is the dominant term, the asymmetry related to the 
A 2 coefficient can be clearly measured, again with a preference for the highest q 2 bins. 
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